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o 

£nQ , A method for a microscopic description of A hypernuclei is formulated in the framework of 

the unitary-model-operator approach. A unitarily transformed hamiltonian is introduced and 
' given in a cluster expansion form. The structure of three-body-cluster terms are discussed 

especially on the A single-particle energy. The A single-particle energies including the three- 
body-cluster contributions are calculated for the US1/2, OP3/2 and Opi/2 states in ]\0, and for 
the 0S]/2, 0p 3/ / 2 , 0p 1; /2, 0d 5 / 2 , 0d 3/ / 2 an d IS1/2 states in \ Ca, using the Nijmegen soft-core 
(NSC), NSC97a-f, the Jiilich A (JA) and JB hyperon-nucleon interactions. It is indicated 
that the three-body-cluster terms bring about sizable effects in the magnitudes of the A 
single-particle energies, but hardly affect the A spin-orbit splittings. 

o ' 

O ' §!• Introduction 

m 

I Much attention has been paid to the study of hypernuclear physics in recent 

years. One of the fundamental problems in hypernuclear physics is to understand 
properties of hypernuclei starting from the basic interaction between two kinds of 
constituent particles, nucleon (iV) and hyperon (Y). Our knowledge on the YN 
interaction is, however, quite inadequate because of the very limited experimental 
information on the YN scattering. In such a situation many-body theoretical stud- 
ies of the hypernuclear structure could provide alternative information on the YN 
interaction. For this purpose, it would be desirable to have a theoretical method 
^ \ which enables us to make a high-precision calculation. 

In many of microscopic structure calculations of hypernuclei, the G-matrix has 
been introduced as a YN effective interaction. The G-matrix approach has been 
widely applied, and the use of this type of effective interaction has made important 
contributions to the description of hypernuclei. ^ " 7 ) However, we may say from a 
formal point of view that the G-matrix is energy-dependent and non-hermitian, and 
does not have the property of decoupling between a model space of low-lying two- 
particle states and its complement. The G-matrix itself is not considered to be an 
effective interaction in a formal sense of the effective interaction theory because it 
does not have the property of decoupling. 8 ) In order to derive such an effective in- 
teraction we should add some higher-order corrections such as folded diagrams. 9 )' 10 ) 
It would be desirable if we could have a theory to describe many-body systems in 



E-mail address: fuji2scp@mbox.nc.kyushu-u.ac.jp 
E-mail address: okamoto@mns.kyutech.ac.jp 
E-mail address: suzuki@mns.kyutech.ac.jp 



typeset using VTpTfX.sty <ver.l.O> 



2 



S. Fujii, R. Okamoto and K. Suzuki 



terms of an energy-independent and hermitian effective interaction with the property 
of decoupling. 

We have been proposed a many-body theory, the unitary-model-operator ap- 
proach (UMOA), u ) that was formulated on the basis of such an effective interaction. 
The UMOA was applied to finite nuclei 16 O n ) and 40 Ca, 12 ) and some encouraging 
results were obtained. Furthermore, in the previous work, 13 ) we extended the for- 
mulation of the UMOA to a description of A hypernuclei and made calculations of 
properties of ^ 7 which is a typical closed-shell nucleus plus one A system, using 
some realistic YN interactions. 

The formulation of the UMOA is based essentially on a cluster expansion of a 
unitarily transformed hamiltonian in terms of a two-body correlation operator. 14 )' 15 ) 
In the previous work we gave a method of treating correlations within the two-body- 
cluster term, and the three-or-more-body-cluster terms remain to be evaluated. It 
has been known that the expansion of the unitarily transformed hamiltonian does 
not, in principle, terminate in finite order. Therefore, it needs to give a method 
of evaluating the three-or-more-body-cluster terms and confirm numerically conver- 
gence of the cluster expansion. It is of high interest to investigate to what extent 
the three-body-cluster term would affect the A single-particle energy, especially, the 
A spin-orbit splitting. The calculation of the A spin-orbit splitting is one of the 
important tasks in hypernuclear physics because it has a possibility of providing 
information to the determination of some parameters in the YN interaction such as 
the strengths of the spin-orbit and antisymmetric spin-orbit components. 

The purposes of this work are to formulate a method of evaluating the three- 
body-cluster effect in the framework of the UMOA, and apply it to and ^Ca 
by employing various realistic YN interactions. This paper is organized as follows: 
In §2 we present a formulation of the UMOA to treat the three-body-cluster term. 
In §3 we apply the UMOA to ^ 7 and 41 Ca by employing the Nijmegen soft-core 
(NSC), 16 ) NSC97a-f, 1? ) the Julich A (JA) and JB 18 ) YN potentials. In §4 we 
make some concluding remarks. 



§2. Formulation 



2.1. Cluster expansion of the unitarily transformed hamiltonian 

We first consider a hamiltonian of a hypernuclear system consisting of nucleons 
and one A (or S) in the second-quantization form as 

H = ^("1^1 1/3)4^/3 + 4 (a/3\v NlN2 \j5)cicy s c 1 

+ + Am Y \v)d\d v + Y (^VYNAv&d^cpdv, (2-1) 

Y=A,E f = A,E 

where c f (c) is the creation (annihilation) operator for a nucleon in the usual nota- 
tion, and cf (d) is the creation (annihilation) operator for a hyperon, A or S. The 
kinetic energies of a nucleon and a hyperon are denoted by t^ and ty, respectively. 
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The vn 1 n 2 an d vyni represent the nucleon-nucleon (NN) and hyperon-nucleon ( YN) 
interactions, respectively. The notations a, (3, 7 and 5 are used for the sets of quan- 
tum numbers of nucleon states, and and v for those of hyperon states. We note 
here that \af3) is the antisymmetrized and normalized two-body NN state. The 
Amy denotes the difference between the rest masses of A and S. 

In nuclear many-body problems, it is important to treat properly the short-range 
two-body correlation because realistic NN and YN interactions have strongly repul- 
sive cores at short distance. For this purpose, we introduce a unitary transformation 
of the hamiltonian as 



where S is the sum of anti-hermitian two-body operators for NN and YN systems 
defined as 



H = e~ s He- 



,5 



(2-2) 



S = S (NN) + S (YN) 



(2-3) 



with 



s (nn) = -J2 {af3\S NlN2 \-f5)clcy 5 Cy, 



(2-4) 




(2-5) 



fiau/3 



According to Providencia and Shakin, 14 ) we make cluster expansion of the 
unitarily transformed hamiltonian as 



(2-6) 



where the first three terms are written explicitly as 




(2-7) 



fJLV 

Y=A,S 





(2-8) 



and 
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Y=A,E 

~ {^\uy Nl \v(3)d\^ a ci3d v . ( 2 ' 9 ) 

liau/3 
Y=A,S 

Here \ct(3~f) is the antisymmetrized and normalized three-body NNN state, and the 
NN part of \(J,a(3) is antisymmetrized and normalized. In the above, we have used 
the notations given as follows: Since the exponent 5 is a two-body operator, the 
one-body operators for k = N\ and Y in H are unchanged and given by 

fiNi = *iVi + u Nl , (2-10) 



hy = t y + uy + Amy ■ (2-11) 

The terms uu for k = N\ and Y are the auxiliary single-particle potentials of a 
nucleon and a hyperon, respectively, and in this stage they are arbitrary The Vki for 
{kl} = {N1N2} and {YNi} are the transformed two-body interactions for NN and 
YN systems, and they are given by 

v Nl N 2 = e~ SN ^ (h Nl + h N2 + v NlN2 )e SN ^ - (h Nl + h N2 ), (2-12) 

vy Nl = e- s ™i (h Y + h Nl + v YNl )e s ™i -(hy + h Nl ). (2-13) 

The transformed hamiltonian H contains, in general, three-or-more-body interac- 
tions even if the starting hamiltonian H in Eq. (2T) does not include three-or- 
more-body interactions. The transformed three-body interactions v^im for {Mm} = 
{N1N2N3} and {YN1N2} are given, respectively, by 

v Ni n 2 n 3 = e~ SN ^ N 3 (h Nl + h N2 + h Ns + v NlN2 + v N2Ns + v NsNl )e SN ^ N 3 

— {hN! + hft 2 + /lAf 3 + VN!N 2 + VN2N3 + VN 3 Nt), (2T4) 

vy Nl N 2 = e~ SYN ^ (hy + h Nl + h N , 2 + v YNl + v NlNa + v N2 y)e SYN ^ 

-(hy + h Nl + h N2 + v YNl + v NlN2 + i)Ar 2 y), (2-15) 

where the operators Ski m for {kirn} = {N1N2N3} and {YN1N2} denote the sum of 
the two-body correlation operators defined as 

Sn 1 n 2 n 3 = 5ViiV2 + Sn 2 n 3 + Sn 3 Ni, (2-16) 



Syn!N 2 = Syn 1 + Sn!N 2 + Sn 2 y- 



(2-17) 
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We here note that the unitary transformation of the auxiliary potentials ut for 
k = N\ (N2) and Y generates two-body operators written as 

u Nl N 2 = e- SN i N 2 (u Nl + u N2 )e SN ^ - (u Nl + u N2 ), (2-18) 

uyn, = e- SYN i(u Y +u Nl )e SYN i - (u Y +u Nl ). (2-19) 

In nuclear many-body problems, it is important how to choose the auxiliary 
potential Uk- In general, it is useful to introduce un 1 (un 2 ) an d U Y as self-consistent 
potentials defined with the transformed two-body interactions vn 1 n 2 an d vyn ± as 

(a\u Nl \P)= J2H\v NlN2 \(30, (2-20) 
(jji\uyW)= E^NJ^). (2' 21 ) 

£<PF 

where pf is the uppermost occupied level, and the notation £ means an occupied 
state for nucleons. If and uy are defined as given in Eqs. (2-20) and (2-21), one 
can prove that the second and fourth terms on the right-hand side of Eq. (2-8) are 
cancelled by the contributions of the bubble diagrams of the first and third terms, 
respectively. In the same way, the second and fourth terms on the right-hand side 
of Eq. (2-9) are also cancelled by the bubble-diagram contributions of the first and 
third terms, respectively. Therefore, the cluster terms and .ff( 3 ) include 

only the one-, two- and three-body operators, respectively, if we write them in the 
normal-product form with respect to particles and holes as discussed in Ref. 11). As 
has been well known, this cancellation mechanism is advantageous in perturbative 
calculations, because there is no need of evaluating u-insertion diagrams and bubble 
diagrams. 

An important problem in the present approach is how to determine the two-body 
correlation operators Ski for {kl} = {N1N2} and {YNi}. These operators are given 
as solutions to equations of decoupling as 

Qkie- Skl (h k + h l + v M )e Skl P k i = 0, (2-22) 

where P^i and Qki for {kl} = {N1N2} and {YN±} are projection operators which act 
in the space of two-body states and project a two-body state onto the low- momentum 
model space and the high-momentum space, respectively. It has been proved that 
Eq. (2-22) for Sm can be solved in a nonperturbative way 11 )' 13 ) under the restrictive 
conditions 

PmSmPm = QkiSkiQki = 0. (2-23) 

In the previous work 13 ) the structure of the two-body-cluster term for the A 
single-particle energy was studied, and three-or-more-body cluster terms were not 
considered. After the evaluation of the two-body-cluster term, the next important 
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term will be the three-body-cluster (TBC) term. The significance of studying the 
TBC term may be understood from the following two reasons: The first is to examine 
whether the TBC effect is large or not in comparison with the one- and two-body- 
cluster effects, in other words, whether the cluster expansion of H is convergent so 
fast or not as to be useful in a practical application. The second is that, when we 
investigate effects of the real three-body force, we must be careful of the possibility 
that such three-body-force effects may interfere with effects of the effective three- 
body force. 

2.2. Structure of the three-body- cluster term 

We here discuss general structure of the TBC term in Eq. (2-9). The TBC 
term contains the transformed three-body NNN and YAW interactions. These 
interactions include the bare two-body AW and YN interactions which have strongly 
repulsive cores. Therefore, we represent the transformed three-body interactions in 
terms of the transformed two-body AW and YN interactions va^a^ an d fyjvi which 
are well-behaved in the actual calculation. 

Using a factorization formula of Campbel and Hausdorf, we have for exp(Sijk) 

e Sijk — e Sij e S(ij)k — e Sjk e S(jk)i — e Ski e S(ki)j (2-24) 

with 

S(ij)k = Sjk + Ski — 2 [Siji Sjk + Ski] 

~ [[^iji Sjk + Ski], Sijk + Sij] + • • -. (2-25) 

Hereafter, a set {ijk} denotes that of {N1N2N3} or {YW1W2}. The correlation 
operators Suk)i an d S(ki)j are a ^ so written in a similar way. Substituting exp(Sijk) 
in Eq. (2-24) into Eqs. (2-14) and (2-15), we have 

= + (2-26) 

with 

4jl = E {e- s W)"Vije s W)" - Vij) (2-27) 

(ijk) 

and 

4jl = E [e- S ^ k (hi + hj)e s W - (h % + hj)] 

(ijk) 

+(hi + hj + h k ) - e~ s ^ k (hi + hj + h k )e s '^ , (2-28) 

where J2(ijk) denotes the summation of exchange terms defined for an arbitrary 
function f(ijk) as 



(ijk) 



(2-29) 
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The terms v^ k and v^l can be expanded into series in powers of Sij, Sjk and 
Ski- The term v^l is given by 

d ijl = 2 S Jk + S ki] ~ ^[Vij, [Sij, S jk + S ki ]} 

(ijk) 

+\[[vij,Sj k + S ki ], Sjk + Ski]} H • (2-30) 

Similarly, v^ k is given by 

d ijl = S i~2^ hi + h j^ s ij],S jk + 

+^[[[^ + hj, S^], Sjk + Sij*]} + ■ ■ •. (2-31) 

The expressions of v^j k and give us desirable forms for the calculation of 
the TBC terms because they do not include the bare interaction Vij and are written 
in terms of only the terms hi, Vij and S^. The TBC effect on the ground-state 
and one-body energies for ordinary nuclei such as 16 and 40 Ca has already been 
studied. n )' 12 )' In the present work we study the effect of the TBC term on the 
energy of A. 

2.3. Derivation of the three-body- cluster contribution to the A single-particle energy 

The contribution of the TBC term to a A single-particle energy in a state \a 
is derived from Eq. (2-9) by making the summation of diagonal occupied states of 
nucleons [In and vn as 

AE^™ C) = ^ (XaHnvn\van 1 n 2 \^aIi-ni< / n) 

- X] {^aVn\uanMmin), (2-32) 



where van^n? ls the transformed three-body ANN interaction Vijk for {ijk} = 
{AN1N2} in Eq. (2-26), and uani is the transformed operator of the one-body po- 
tential defined for Y = A in Eq. (2-19). 

In order to obtain the TBC effect in Eq. (2-32), we need to calculate the matrix 
element of the two-body correlation operator Sij, the two-body effective interaction 
Vij and the single-particle potential The calculation procedure for these quantities 
has already been given in Refs. 11) and 13) in detail. Therefore, we here present 
only the outline of the calculation procedure. In an actual calculation for the YN 
system, we employ a two-step procedure by introducing two kinds of model spaces, 
namely, large and small model spaces as shown in Figs. 1 and 2. We here use the 
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harmonic oscillator wave functions as basis states. The numbers pp, px, Pi and P2 
in Figs. 1 and 2 are defined as 

pp = 2n N + l N , (2-33) 

P2 = 2n A + l A , (2-34) 

px = 2n Y + W + 2n N + l N , (Y = A, S) (2-35) 

and 

pi = 2ua + Ia + 2n N + l N , (2-36) 

where n& and Ik for k = A, S and N are the harmonic oscillator quantum numbers. 
The number p2 is introduced so as to specify the uppermost bound state of A. 

In the first-step calculation, we calculate the AN effective interaction v^ j^ acting 
in the large model space P^h by solving the equation of decoupling between the P^ 
and the - Q^v) + (Qen ~ ^S)} s P aces in E Q- O 22 )- 13) Jt is noted here 

(X) (X) 

that the Q A ^ and Q SI ^ spaces are the excluded spaces due to the Pauli principle 
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(a) 

Fig. 3. Calculated diagrams in the second order for the A single-particle energies. The AN effective 
interaction van is represented by the wiggly-line vertex, and the vertex x with the dashed line 
represents the non-diagonal part of the one-body hamiltonian hk = t^+Uk for k = A and N. The 
propagations of N and A are represented by the single and double external lines, respectively. 

for nucleons. The A potential energy in the first step is calculated as given in 
Eq. (2-21), using the AN effective interaction 

e (i) 



AN- 



(2) 

In the second-step calculation, the AN effective interaction v AN acting in the 

(2) (2) 

P\x space and the correlation operator S A ^ are obtained by solving the equation 



of decoupling between the P A ^ and the Q AN spaces, using v A ^ and u K A . The A 



(2) 



„~,(1) 



(2) 

potential energy u A is finally calculated in Eq. (2-21), using v A ^. However, at this 

(2) (2) 

stage, we consider only the diagonal part of h A = tj[ + u A as the unperturbed 
one-body hamiltonian. Therefore, we evaluate perturbation corrections arising from 

(2) 

the non-diagonal terms of h A as shown in diagram (a) of Fig. 3. Furthermore, the 



(2) 



effective interaction between AN states in the P A ^ and P A ^ spaces induces core 
polarization in the closed-shell nucleus. We take account of this effect by calculating 
diagrams (b) and (c) of Fig. 3, using the AN effective interaction v A lf determined 
in the first step. The validity of this two-step calculation procedure is confirmed 
numerically in the previous work. 

In order to calculate the TBC effect on the A single-particle energy, we need 
the effective interaction and the correlation operator for the NN system which are 
determined by solving the equation of decoupling in Eq. (2-22). Using the two-body 
effective interactions and the two-body correlation operators for the AN and NN 
systems, one can derive that an approximate expression of the TBC terms from 
Eq. (2-32) as 



,00 



AE 



a™ '* — 9 {^a^n^nI^ H^Wl, Sn 1 n 2 + Sn 2 a], Sn 1 n 2 + Sn 2 a]\^aHnvn) 



[l N U N <p F 

(ANlN2) 



(or) + (6t) + Ot) + (<2t) + (e T ) + (/t), 



(2-37) 



where 



(ot) 



PN^N<PF a N>PF 
"A<P2 a A>P2 



a{AN) c(AN) 

^A^N^ApN v APNOtA<*N Ot A OtN^AVN ' 
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(hr) = - V V V (AN) S {NN) S (NN) 

v T / 2 ^ ^APN^APN Pnvn(*nPn ajV/3iVMiV^JV' 

PN^NPN<PF ajv/3jV>PF 

(r^ = -2 V V f; (/1#) <7 {WV) <7 {ylA ° 

V \ A a N a A PN UnvncinPn a A f3 N \ A v N > 

Pnvn<Pf a N fS N >p F 
a A >p2 

(drr) = V V 5 (JVJV) S (AN) S (AN) 

V ^/^ 2_^i PnO-NVnPn ^A^NCtAQ-N OlaPn^AVn' 1 

Pn^n<Pf a N (3 N >p F 

CtA>P2 

(Wl = V V v (AN) S (AN) ? {AN) 

V ^ u a A v N j3 A n N kD \ A ^ N a A a N (3 A a N \ A v N i 

Pn"n<Pf ojv>Pf 

(fcO = -- V V v (AN) ? (NN) S (NN) (2-381 

wT; ^ A/1QJVA/17JV p, N v N a N (i N inPn^nvn" v ' 

Pnvn<Pf ajv/3jv7iV>PF 

In the above, we have used the notations defined by 

vSUnSn = ( a NPN\v Nl N 2 \lNSN) (2-39) 

and 

S Z N plyN6 N = (<xnPn\S Ni n 2 \inSn)- (2-40) 

The contributions of the TBC terms to come from the effective three- 

body interactions v^ k and v^. We here note that the first and second terms of 
v\jl in Eq. (2-30) do not bring about non-zero contributions to the A single-particle 
potential. On the other hand, the contributions of vf^ k in Eq. (2-31) vanish up to 

second order in Sij. Therefore, the contributions of v^jj. and v^j, start with the 
second and third order terms in Sij, respectively. These facts concerning the TBC 
terms can be proved by using the property of decoupling of Vij and the restrictive 
conditions for Sij in Eq. (2-23) as discussed in Ref. 15). On the assumption that 

the matrix elements of Sij are sufficiently small, we neglect the contribution of v^.. 
The validity of this assumption will be confirmed numerically in §3. 
The A single-particle energy is given finally by 

Ex A = (XaMXa) + (\a\u { ?\\ a ) + AE^ A + AE { ™ C \ (2-41) 

where the first term is the kinetic energy, the second the first-order potential energy, 
the third the sum of the perturbation corrections given in Fig. 3 and the last the 
sum of the TBC contributions given in Eq. (2-37). 

2.4. Diagrammatical expression of the three-body-cluster term 

We here make some discussions on the relation between the UMOA and the 
G-matrix theory in the calculation of the A single-particle energy in order to clarify 
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Fig. 4. Diagrammatical expression of the first-order A potential energy in terms of the AN effective 
interaction van obtained in the UMOA and the G-matrix Gan- 



Fig. 5. Diagrams corresponding to the leading order terms of the three-body-cluster terms in terms 
of the G-matrix shown by the wavy line. 

the physical meaning of the TBC effect. The formal relation between the effective 
interaction in the UMOA and the G-matrix has been discussed in Refs. 11) and 15). 
The AN effective interaction van in the UMOA in the calculation of the first-order 
A potential energy ua can be given in an expansion form in terms of the G-matrix 
as shown diagrammatically in Fig. 4. The second diagram on the right-hand side in 
Fig. 4 is a kind of folded diagram. 9 ) The leading contributions of the TBC terms 
in Eq. (2-38) can also be given diagrammatically in terms of the G-matrix as shown 
in Fig. 5. It is noted that the second diagram on the right-hand side in Fig. 4 
has a minus sign which is against the diagram rule. However, this diagram is just 
cancelled by diagram (ar) in Fig. 5. It should be noted that such a folded diagram 
emerges again as one of higher-order corrections with the vertices of the transformed 
interaction v^. 




'an 



+ 
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§3. Numerical results for A O and ^ Ca using the Nijmegen and the 

Jiilich YN potentials 

We performed calculations employing the Nijmegen soft-core (NSC), 16 ) NSC97a- 
f, 17 ) the Jiilich A ( JA) and JB 18 ) potentials for the YN interaction. We used the 
harmonic oscillator basis states and took a common fifi to single-particle states of 
A, £ and N. We used the values %S2 = 14MeV for \ 7 and 12MeV for ^Ca. The 
dependence of the calculated effective interaction and the single-particle energy on 
the value %fi has already been examined, and we confirmed that the dependence 
was quite small around %Q = 14MeV in 16 15 ) if perturbation corrections were in- 
cluded. As for the nucleon single-particle potential un, we used the fixed data 11 )' 12 ) 
calculated using the Paris potential, 19 ) and we have assumed that us = 0. The 
angular momentum I of the partial wave of the NN and YN interactions is taken 
into account up to I = 4. The numbers p2, Pf, Pi and px are taken as p2 = 1, Pf = 1, 
Pi = 8 and px = 12 for \ 7 0, and p2 = 2, pp = 2, pi = 10 and px = 12 for ^Ca. 

The formulae for calculating the TBC contributions to the A single-particle 
energy have been derived in Eq. (2-38). The TBC terms are given in powers of the 
two-body correlation operator Sij as given in Eqs. (2-30) and (2-31). In the UMOA, 
it is necessary that the matrix element of the correlation operator is sufficiently 
small in the sense of cluster expansion. In Table I typical matrix elements of Sij 
for the NN and AN systems are given. We present only the matrix elements of 
S^ of comparatively large magnitudes among all the matrix elements of Sij. The 
magnitudes of the matrix elements are at most of the order of 10~ 2 . Besides, the TBC 
effect yields the contributions of the order of the square of Sij as shown in Eq. (2-38). 
Thus, we expect that the TBC contribution is considerably small compared with the 
contribution from the one- and two-body-cluster terms. 

We calculated the A single-particle energies with the correction terms (a), (b) 
and (c) given in Fig. 3 and the TBC terms (ar)-(fT) i n Eq. (2-38) for l j0 using 
the NSC97f potential as shown in Table II.*) One see that the TBC contributions 

*-* The mark *** in Tables II and III means that the calculated absolute value is less than 
O.OlMeV. 

Table I. Typical matrix elements {aa\S^ N \bb) j,t and (aa\S^\bb) ,/,t=i/2 hi ^ 7 for the two-body 
NN and AN correlation operators. The Paris and the NSC97f potentials are used as the NN 
and the YN interactions, respectively. The letters a and b denote single-particle orbits. The J 
and T mean the total angular momentum and isospin, respectively, for two-body states. 



a 


b 


J 


T 


(aa\Sxh\bb)j lT 


T 


{aa\S ( ^\bb)j,T 


lSl/2 


0Si/2 


1 





2.32 x 10~' A 


1/2 


5.38 x 10" 3 


1*1/2 


0Si/2 





1 


1.69 x 1(T 2 


1/2 


9.42 x 1(T 3 


1*1/2 


Opi/2 


1 





-2.43 x 1(T 2 


1/2 


2.49 x 1(T 3 


1*1/2 


Opi/2 





1 


-5.75 x 1(T 3 


1/2 


9.36 x 1(T 4 


0rf 3 /2 


OSl/2 


1 





8.31 x 1(T 3 


1/2 


1.51 x 1(T 4 


0d 3 / 2 


0si/2 





1 


7.82 x 1(T 3 


1/2 


4.58 x 1(T 3 


0rf 3 /2 


0pi/2 


1 





-6.48 x 1(T 3 


1/2 


3.67 x 1(T 3 


0d 3/2 


0pi/2 





1 


-8.85 x 10" 3 


1/2 


1.29 x 1(T 3 
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are considerably small in comparison with the contributions from the one- and two- 
body-cluster terms. This ensures that the cluster expansion method in powers of the 
correlation operator Sij may be justified in the calculation of the A single-particle 
energy. 

One also see that two of the TBC terms, (bp) and (fx), have dominant con- 
tributions among the six TBC terms. It may be clear from Fig. 5 that these two 
TBC terms include the vertices inducing typical two-particle two-hole excitations of 
the core nucleus in intermediate states. On the other hand, the TBC terms (ax) 
and (cx)-(ex) do not include this type of excitations. The two terms (bx) and (fx), 
however, give the contributions of the opposite sign. Therefore, the sum of the TBC 
terms becomes rather small. The effect of the TBC terms is, on the whole, repul- 
sive and at most about 3-4% of the A one-body potential energies of the one- and 
two-body-cluster terms for the Op and Os states in ^ O as given in Table II. 

The calculated value of the spin-orbit splitting of the A single-particle energies 
for the Op orbits in ^ 7 is l.llMeV for the first-order potential energy. Corrections 
(a), (b) and (c) in Fig. 3 yield the contribution of — O.llMeV to the A spin-orbit 
splitting. The sum of the TBC terms for the A spin-orbit splitting is only — 0.06MeV. 
We may conclude from the present calculation that the magnitude of the A spin-orbit 
splitting is determined dominantly by the contribution of the one- and two-body- 
cluster terms, and the polarization of the core nucleus does not play an important 
role. This trend is also confirmed for the other YN interactions used in this study. 

The calculated result in ^Ca for the NSC97f is shown in Table III. The effects 
of the second-order and the TBC correction terms on the A single-particle energy 
are not so different from those in \ 7 0. The magnitude of the TBC terms is at most 
3% of the A one-body potential energy of the one- and two-body-cluster terms for 
the Os state in ^Ca. We note here that the A spin-orbit splittings are calculated for 

Table II. Contributions to the A single-particle energies in ^ 7 for the NSC97f potential. KE and 
PE stand for the kinetic and first-order potential energies, respectively. Rows (a), (b) and (c) 
are the contributions of diagrams (a), (b) and (c) in Fig. 3, respectively. Rows (aT)-(fr) mean 
the contributions given in Eq. (2-38). All entries are in MeV. 





Osi/2 


0P3/2 


Opi/2 


Aei s 


KE 


10.50 


17.50 


17.50 




PE 


-23.13 


-17.20 


-16.09 


1.11 


(a) 


-1.84 


-2.35 


-2.42 


-0.08 


(b) 


-1.20 


-0.37 


-0.13 


0.24 


(c) 


-0.55 


-0.98 


-1.25 


-0.28 


Subtotal 


-16.22 


-3.40 


-2.40 


1.00 


(a T ) 


-0.01 


-0.01 


-0.01 


*** 


(b T ) 


1.49 


1.22 


1.15 


-0.07 


(ct) 


-0.09 


-0.16 


-0.16 


*** 


(dx) 


*** 


*** 


-0.02 


-0.02 


(ex) 


0.01 


0.02 


0.02 




(fx) 


-0.43 


-0.46 


-0.44 


0.02 


Sum of TBC 


0.97 


0.60 


0.54 


-0.06 


Total 


-15.25 


-2.80 


-1.86 


0.95 
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both the Op and Od states in 4 Ca. The values of the A spin-orbit splittings for the 
Op and Od states are 0.79MeV and LllMeV, respectively. The net contribution of 
core polarizations to the A spin-orbit splittings is very small for both the Op and Od 
states. This situation is quite different from that in the ordinary nucleus 40 Ca, 12 ) 
as well as in 16 0. n ) It has been observed that many-body correlations bring about 
corrections of more than 30% for the spin-orbit splittings of nucleon states. 

It might be of interest to compare the results obtained in this study with those 
in other many-body methods. Recently, the A spin-orbit splittings in were 
calculated on another method by Motoba 20 ) in which configuration mixing of the 
A single-particle states is taken into account, using the YNG interaction 21 ) derived 
from the NSC97 potentials. The result for the NSC97f was almost lMeV for the Op 
states. This value is consistent with our result of 0.95MeV for the NSC97f as given 
in Table II. 

We now discuss the results obtained for various YN interactions. The partial- 
wave contributions of the AN effective interactions to the first-order A potential 
energy are given for various YN potentials, namely, the JA, JB, the NSC*) and 

*-* The mark (*) for the NSC in Table V and Fig. 6 means that the antisymmetric spin-orbit 
force is not included. This potential was also used in the calculation of properties in in the 

Table III. Results of fCa for the NSC97f. The notation is the same as in Table II. All entries are 
in MeV. 





0si/2 


0P3/2 


0pi/2 


0d 5 /2 


0^3/2 


lSl/2 


KE 


9.00 


15.00 


15.00 


21.00 


21.00 


21.00 


PE 


-28.01 


-22.75 


-21.76 


-18.08 


-16.63 


-16.36 


(a) 


-2.16 


-3.19 


-3.25 


-3.69 


-3.83 


-2.96 


(b) 


-1.44 


-0.99 


-0.87 


-0.39 


-0.15 


0.16 


(c) 


-0.58 


-1.02 


-1.25 


-1.32 


-1.73 


-1.62 


Subtotal 


-23.19 


-12.96 


-12.14 


-2.48 


-1.33 


0.22 


(ax) 


*** 


*** 


*** 


*** 


*** 


-0.01 


(b T ) 


1.46 


1.30 


1.24 


1.12 


1.04 


1.03 


(ct) 


-0.05 


-0.08 


-0.08 


-0.13 


-0.12 


-0.15 


(d T ) 


*** 


*** 


*** 


*** 


0.02 


0.01 


(ex) 


0.01 


0.01 


0.01 


0.02 


0.02 


0.02 


(fx) 


-0.38 


-0.44 


-0.44 


-0.48 


-0.46 


-0.45 


Sum of TBC 


1.03 


0.78 


0.74 


0.53 


0.49 


0.46 


Total 


-22.16 


-12.18 


-11.39 


-1.95 


-0.84 


0.68 



Table IV. Partial-wave contributions to the first-order A potential energy for the Osi/2 state in 
"O. All entries are in MeV. 



I : even I : odd 





^0 


3 ft 




x Pi 


'Pi 


J P 2 


Total 


NSC97a 


-2.37 


-22.93 


-0.03 


0.84 


0.80 


-0.98 


-24.68 


NSC97b 


-3.52 


-22.35 


0.04 


0.90 


0.98 


-0.88 


-24.82 


NSC97c 


-5.06 


-22.00 


0.18 


0.95 


1.12 


-0.78 


-25.60 


NSC97d 


-7.16 


-20.55 


0.30 


1.06 


1.40 


-0.56 


-25.52 


NSC97e 


-8.42 


-19.28 


0.35 


1.14 


1.62 


-0.41 


-24.99 


NSC97f 


-9.46 


-17.11 


0.32 


1.29 


2.00 


-0.14 


-23.10 
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NSC97a-f potentials for the Osi/2 state, in Table IV for ^ O, and in Table V for 
^Ca. It is seen that the partial- wave contributions are very different, especially in 
the spin singlet and triplet states, dependently on the YN potentials employed. It 
is noted, however, that the ratios of the contributions of the spin singlet and triplet 
states vary smoothly from the NSC97a to the NSC97f. This variation is caused 
mainly by the difference in the strength of the spin-spin interaction of the NSC97 
potentials. The similar tendency was observed in the calculation of the potential 
energy of A contained in nuclear matter. 17 ) 

One of the features of the NSC97 model for the YN interaction is that the 
NSC97a-f potentials are given by varying smoothly the magnetic Fj (F + D) ratio 
a™ for the vector mesons in the NSC97a-f potentials. This change in the magnetic 
ratio yields the variation of the strength of not only the spin-spin interaction but also 
other components including the spin-orbit interaction. According to the variation of 
the YN interactions, the calculated result of the A single-particle potential energies 
vary, dependently on the YN potentials employed as seen in Tables IV and V. It is 
important to investigate the difference among the NSC97 potentials quantitatively in 
the calculation of properties in finite A hypernuclei. This quantitative study would 
help to improve the YN interaction, complementarily to the progress in relevant 

previous work. 13 - ) On the other hand, the NSC newly used in the present study is the complete 
version of the NSC potential code which includes the antisymmetric spin-orbit force. 

Table V. Partial- wave contributions to the first-order A potential energy for the Osi/2 state in * x Ca. 
All entries are in MeV. 









I : 


even 












i D 1 


'D 2 


*D 2 




JA 


-3.55 


-23.76 


-0.12 


-0.13 


-0.14 


-0.15 


JB 


0.10 


-31.71 


-0.09 


-0.13 


-0.14 


-0.11 


NSC97a 


-2.80 


-29.42 


-0.09 


-0.17 


-0.22 


-0.21 


NSC97b 


-4.40 


-28.69 


-0.08 


-0.19 


-0.23 


-0.22 


NSC97c 


-6.57 


-28.25 


-0.08 


-0.20 


-0.23 


-0.23 


NSC97d 


-9.55 


-26.27 


-0.08 


-0.23 


-0.24 


-0.24 


NSC97e 


-11.33 


-24.52 


-0.07 


-0.24 


-0.25 


-0.24 


NSC97f 


-12.82 


-21.55 


-0.05 


-0.26 


-0.27 


-0.24 


NSC 


-13.83 


-8.92 


-0.06 


-0.25 


-0.23 


-0.19 


NSC(*) 


-13.12 


-9.50 


-0.05 


-0.24 


-0.23 


-0.19 






I : 


odd 










"Po 


'P 


d P 


6 P2 




Total 


JA 


0.76 


0.82 


1.41 


-0.31 




-25.17 


JB 


0.70 


1.16 


2.00 


0.52 




-27.71 


NSC97a 


-0.05 


1.70 


1.75 


-2.02 




-31.53 


NSC97b 


0.09 


1.83 


2.14 


-1.83 




-31.57 


NSC97c 


0.38 


1.93 


2.42 


-1.63 




-32.48 


NSC97d 


0.63 


2.16 


3.02 


-1.21 




-32.01 


NSC97e 


0.72 


2.34 


3.47 


-0.90 




-31.03 


NSC97f 


0.67 


2.67 


4.26 


-0.36 




-27.97 


NSC 


0.49 


1.79 


2.19 


-3.67 




-22.69 


NSC(*) 


0.51 


1.87 


2.25 


-3.71 




-22.42 
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Fig. 6. Dependence of the calculated A single-particle energies in ^ 7 on various YN interactions. 
All entries are in MeV. 



experiments of A hypernuclei. 

In Figs. 6 and 7 we show the calculated A single-particle energies including the 
second-order perturbation corrections and the TBC terms for the JA, JB, the NSC 
and NSC97 interactions. It is noted here that the calculated results for the NSC97f 
potential are, on the whole, fairly good in comparison with the empirical values for 
^ 6 and A °C&. 7 ) Our result for the NSC97f is consistent with the fact that the 
calculations of few-body A hypernuclei such as ^H 17 ** and ^He 22 ) using the NSC97f 
can reproduce successfully experimental energy levels. 

The calculated A spin-orbit splittings of the Op states vary in the range of 0.49- 
1.31MeV from the NSC97a to the NSC in x j0 as shown in Fig. 6. Similarly, the A 
spin-orbit splittings of the Op and Od states vary in the range of 0.35-1.08MeV and 
0.61-1.55MeV, respectively, in ^ Ca as shown in Fig. 7. The experimental values 
of the spin-orbit splitting of a nucleon are known as 6.20MeV for the Op states in 
16 and 6.75MeV for the Od states in 40 Ca. Therefore, the obtained values of the A 
spin-orbit splitting are considerably smaller than those of a nucleon. Recently, some 
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Fig. 7. Dependence of the calculated A single-particle energies in ^Ca on various YN interactions. 
All entries are in MeV. 



experiments in connection with the A spin-orbit splittings have been performed at 
KEK for fY and fV 23 ) and at BNL for fC 24 ^ and ^Be. 25 ) Higher resolution 
has been attained in these experiments than that in old experiments though the 
analysis of these new experiments are still being in progress. Therefore, it will 
be highly desirable that the magnitude of the A spin-orbit splitting is established 
experimentally. 

§4. Concluding remarks 

A method of evaluating the three-body-cluster (TBC) effect was presented for 
the calculation of A single-particle energies. The TBC terms were generated as 
the transformed three-body interactions among the YNN and NNN systems. The 
transformed three-body interactions were given in the cluster expansion form in 
powers of two-body correlation operators of the YN and AW systems. Therefore, it 
is necessary for fast convergence of the cluster expansion that the matrix elements of 
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the correlation operator should be sufficiently small. In our calculations, it was made 
sure that the matrix elements of the correlation operator were quite small, and thus 
the TBC contributions to the A single-particle energies were considerably small in 
comparison with the contributions from the one- and two-body-cluster terms. This 
ensured that the cluster expansion method in terms of the correlation operator was 
acceptable in the calculation of A single-particle energies. 

The present results showed that A in the 0p 1 / 2 state in ^ 7 was not bound for 
the JA and the NSC potentials but bound for the JB and the NSC97a-f potentials. 
The spin-orbit splittings of the A single-particle orbits in the Op states, including the 
TBC contributions, were given by 1.35, 1.49, 1.31 and 0.49-0.95MeV for the JA, JB, 
the NSC and NSC97a-f potentials, respectively. The value of 1.49MeV for the JB 
potential is very large compared with that given in Ref. 7). On the other hand, the 
value 0.95MeV for the NSC97f potential is consistent with that given in Ref. 20). The 
calculated A spin-orbit splittings depend largely on the YN interactions employed. 
The dependence of the calculated results on the YN interactions was also examined 
for ^Ca. It seemed that the NSC97f potential yielded better results for ^ 7 and 
y^Ca in comparison with the relevant experimental values. The present results were 
consistent with those obtained in the calculation of few-body A hypernuclei. 

In the application of the UMOA to A hypernuclei, some problems still remain. 
First, in the present calculation, single-particle energies of E in intermediate states 
were neglected. However, the single-particle spectrum of E has a possibility of 
bringing about some contributions to the AN effective interaction. Although single- 
particle energies of E has not yet been known well, it would be important to inves- 
tigate how much the spectrum of E affects the determination of the AN effective 
interaction. Second, there would be some relativistic effect for the A single-particle 
energy as well as in calculations of ordinary nuclei. In the present approach, the 
relativistic effect has not yet been considered in studies of ordinary nuclei and hy- 
pernuclei. On the other hand, the relativistic many-body approach such as the 
Dirac-Brueckner-Hartree-Fock approach 26 ) does not fully take account of the effect 
of three-or-more-body correlations. Therefore, it would be of great importance for 
the understanding of nuclear properties from the many-body theoretical point of 
view that the relativistic effect is incorporated in the framework of the UMOA. 
This problem, which would be closely related to the treatment of the three-body 
force, 27 )' 28 ) is also a remaining task in the UMOA. 
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